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I. INTRODUCTION 

The light front Tamm-Dancoff (LFTD) method |]J has been introduced as an alternative 
tool to lattice gauge theory to investigate relativistic bound states nonperturbatively. The 
LFTD is a Tamm-Dancoff (T-D) approximation || applied to field theories quantized on the 
light front coordinate. In the usual perturbative field theory quantized in the equal time 
frame, a vacuum state is an infinite sea of constituents such as electrons and photons in QED 
and quarks and gluons in QCD. Bound states of a hydrogen atom in QED, and mesons and 
baryons in QCD, arise as excitations of this sea. 

The T-D method independently considered the possibility of describing the vacuum and 
the bound states with a finite number of particles, and solving a set of coupled integral 
equations. This method was applied to a variety of problems in strong interactions, but it was 
unsuccessful because a large number of amplitudes was required to solve a given problem. 

On the other hand, in the LFTD method all constituents have non-negative longitudinal 
momenta defined by p + = p ^ in the light front coordinate. Then the vacuum of the system 
under consideration can not have constituents, or the physical vacuum is equivalent to the 
bare vacuum in the light front coordinate. Therefore, we may expect this approach to remove 
a serious problem in the T-D approximation in the equal time frame. 

It was found in several problems such as the hydrogen atom, positronium, and the two- 
dimensional Yukawa interaction, that one needed an additional T-D amplitude to obtain the 
particle spectrum ||. We have provided an argument |3| for why this happens to be the case. 

In 1974, 't Hooft || introduced a model, 2 dimensional QCD, that continues to be studied. 
We examine this model in the framework of the LFTD method. 
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II. MESONS IN 2D QCD 

We consider the quantized version of Einstein's equation 

2P(H + H I )\^f >= M 2 |^ >, (2.1) 

in two dimensional QCD with SU(N C ) symmetry, for states with mass M on the light front 
coordinate. Here P is the momentum operator in the SU(N C ) gauge, Hq is the free part of the 
Hamiltonian, given by H = Pj ree + P s ^/, and Hi = P ~ + P 2 is the interaction Hamiltonian. 
We consider a meson state |^(/, f',p) > with momentum p, which is described by 

!*(/,/'; p) >= 

I /S= ~^= S (P ~ g 1 ~ ga)^2(/> gi» 92)& f (/> gi> c^Cf', qa, c)|0 > 

r°° d qi dq 2 dq 3 dq 4 . 
+ 2^ / /o ~7n=- /o — -^=Q(p-gi-g2- g3~g4j 

fl,hJ 3 ,f4Cl,^3,C4 J ° V^h V27Tg2 V2^g3 V2^4 
X ^(Z, /'; /l, /2, /a, /4! gl, g2, g3, g4! Ci, C 2 , C 3 , C 4 ) 

x 6+ (fx, qi,cxp(f 2 ,q2, cz)d(f 3 , q 3 , c 3 )S (/ 4 , g 4 , c 4 )|0 > . (2.2) 



The first order LFTD approximation results from putting Eq. ( [2.2| ) into Eq. (|2.1| ) and project- 
ing the resultant equation onto a state with a fermion and an anti-fermion, and a state with 
two fermions and two anti-fermions. This leads to 

< fi, gi, ci; / 2 , g 2) Ca| (M 2 - 2P# ) |tf >=< / 1? 9l , Cl ; / 2 , g 2 , c 2 |2PF 7 |^ >, (2.3) 

and 

< /i, gi, ci; / a , g 2 , c 2 ; / 3 , g 3 , c 3 ; / 4 , g 4 , c 4 | (m 2 - 2PPT ) 1^ > 
= < A,gi,ci;/a,g2,c 2 ;/ 3 ,g3,C3;/ 4 ,g 4 ,c 4 |2Piy/|* >, (2.4) 



where 



\f f . ^(/l,gl,Cl) rf t (/ 2 ,92,c 2 ) | , _s 

| /l ,g 1 ,c 1 ; /2 ,g 2 ,c 2 >=^= ^^^|0>, (2.5) 



and 



fi, 9i, ci; / 2 , q 2 , c 2 ; / 3 , 93, c 3 ; U, g 4 , c 4 > 

9l, c i) &t (/2, 92, c 2 ) S(f 3 , 93, C 3 ) ^(/k, 94, c 4 ) 



|0 > . (2.6) 



V27T91 V27T92 ^2^93 V27T94 

It is straightforward to derive coupled integral equations for the amplitudes ip 2 and ■0 4 . 
The results are 

S fif S f2f'S cl c 2 F(M; /, /'; g b q 2 )ip 2 (f, /'; 9i, 92) 
N c 2 - 1 £>9 2 Z" 00 ^ 



r r r ^ c ~ L P9 f°° dli , ^ j- 7, , , \ 

2P3 2 X ^ x - ^ dk 

L TcidMp J. 00 '^ Jo q 



j: flavors d\,d2,ds,:colors 



x ^4(/, /'; [j, A], [j, / 2 ]; [9, 9i + fc], [-9 - fc, 92]; [di, tfe], [ds, c 2 ]) 



Z, Z, T d 3 dr,d 2 c 2 p ~j~2 / "9 



^ j: flavors di,d2,d'i:colors 00 ^ 

x ^ 4 (/, /'; [j, [/ 2 , ]]; [9, 9i], [92 -k,-q + k]; [di, ci], [d 2 , d 3 ]), (2.7) 

where p stands for principal value integral, iV c is the order of the gauge group, and we have 
used 

%d;ef = ~ {ScfSde - — <5 C A/} , (2.8) 

F(M;f 1 J 2 ...;q 1 ,q 2 ...) 

= M 2 -2( qi + q 2 + ...)\ h W> + h m + , (2.9) 

I ^9i ^92 J 

^(Z, /'; LA, / 2 ], [/ 3 , A]; [91, 92], [93, 94]; [ci, c 2 ], [c 3 , c 4 ]) 

= ^4(/, /'; /1, /2, / 3 , A; 9i, 92, 9s, 94; ci, c 2 , c 3 , c 4 ) - (1 <-> 2) - (3 <-> 4) 

+ (1 2, 3 <-> 4), (2.10) 

and 

^(■W; /1, /2, /s, A; 91, 92, 93, 94)^4(7, /'; [A, / 2 ], [js, A]; [91, 92], [93, 94]; [ci, c 2 ], [c 3 , c 4 ]) 



2pg 



2tt 



{ S fif S f2f 3 ^Uf'' I c 1 c 4 ;c 2 c 3 ^2( y f, /'; 9i + 92 + 9s, 94) • t- - ^ 2 
(1 <-> 2) - (3 «-> 4) + (1 <-> 2, 3 <-> 4)} 

{dfifJfzftihf'^caincMf, /'; 92, 9l + 93 + 94) 



(?1 + 94) 1 



- (1 <-> 2) - (3 «-> 4) + (1 <-> 2, 3 ^4)} 



(2.11) 



The a term in Eq. ( |2.9|) denotes the meson self energy and is found to be equal to 
-(N 2 -l)g 2 /(2N c n). 

Hereafter, we restrict ourselves to the one flavor case. We define 

P = qi + q2, qi=(^+x^p, ^2(91,92) = b (x), a=^ + x and b = i - x. (2.12) 



Substituting Eq. ( $J% ) into Eq. (|2?7| ), we have 



M 



2 4(m 2 + a) 



1 - 4a; 2 



6 (x) = I (x) + + I 2 {x) 



(2.13) 



where 



and 



io(s) 
h(x) 



-2Nip f %b {x + y), 
J -a y 2 



f° dy f-y 
/ — / dz 



a y 2 Jo F(M;z,a + y,-z -y,b) 



[a — Zf y" 
n 



(b + z) 2 (l + yf 
n 1 



(a-z) 2 (i + y y 



(2.14) 



- {(6T^ + 7W a; + ^-{7^W + 7TT^}^-^' ^ 



l-<(r) = -M 2 p ^ I 'd: 



y 2 Jo F(M; —z + y,a,b — y, z) 



(6-z) 2 (1-y)* 
// 1 
I2 



n 



2/ - zj 



(2.16) 
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Here we use the shorthand notation 

M = (iV c 2 -l)#7(4iV c 7r), Mi = (N c -l)g 4 /(4N c ir 2 ), and n = (N c +l)/N c . 

(2.17) 

With the definition of Eqs. (|2.15|) and ( p.!6| ) we can show 

I 2 {-x) = ±h{x), (2.18) 

corresponding to the symmetry property of bo(±x) = ±bo(x). It is not possible to solve 
Eq. ( |2.13| ) analytically. One usually selects a wavefunction b (x) and tries to match both sides 
of Eq. (|2.13D. It happens that the results do not depend sensitively on the form of b Q (x). In 



other words, different forms of bo(x) lead to similar predictions. Further, the contributions of 
the integrals h(x) and h(x) are at most 10% of that of Iq(x). Suppose we take the simple form 
of power series expansion for 6 - Since we wish to determine the ground state wavefunction, 
which has no node (and thus no odd powers of x), we write 

b Q (x) = 1+Ax 2 + Bx 4 + Cx 6 + ---. (2.19) 



The form of Eq. ( |2.19p is motivated by our assumption that bo(x) is an even function 



and we normalize the first term to 1. The power series is expressible in terms of a linear 
combination of Airy finctions that are relevant to linear potentials. It is surprising, as we shall 
show, that most of the information is present in the first terms 1 + Ax 2 of ( |2.19 ) and that 



higher order terms such as the B and C terms provide a small correction. We expect A to 
be negative in order that bo(x) decrease away from x — 0. The wavefunction 1 — \A\x 2 is a 
maximum at x = and is zero at x = ±|A| -1 / 2 . 

As a warmup, consider the simplest model || of ( 2.13 ) with only Jo (a:) and 

b (x) = 1 + Ax 2 . (2.20) 



We study this model because it contains the ingredients of the resulting relation of the ground 
state mass M and the quark mass m. The additional contributions of Ii(x), h(x) and higher 
order terms such as B and C terms provide a small correction to this relation. 



Define 



m 



M = M 2 /2Nx m 2 = ( ) - 1 

1 1 v 2ivy 



where the relation a = —2Ni is used, and write the simplified (|2.13|) as 

{(1 - 4x 2 )M 2 - 4m 2 } (1 + Ax 2 ) = 4 + A i^Ax 2 - 2(1 - 4x 2 )xln - 1 + 4x 2 | 

From equating the coefficients of the x° and x 2 terms, we obtain 

M 2 — Am 2 = 4 — A, 



(2.21) 



(2.22) 



(2.23) 



4M 2 - (M 2 - Am 2 ) A = -1QA 



(2.24) 



We find that when 6o(^) = 1> Eqs. ( |2.23| ) and ( ^.24| ) yield M = m = 0, a relation we regard 
as a boundary condition. When we eliminate A from (|2.23|) and (|2.24|), we obtain 



M/g = \/6m/g = 2.A5m/g 



(2.25) 



and the condition M be real is m/g < 2/V37T- Also, one can express M 2 and m in terms of 



-4. 



M 2 3 + A/A 



(2.26) 



m 2 1/2 + A/16' 

We return to Iq(x) and I\(x) given by Eqs. (2.14) and (2.15), respectively. The Io(x) is 
written in the form 



Io(x) 



iA/i 



1 -Ax 2 
+ 25ATi| 

+ 20M 



2A/Vi 



Ax' 



I -Ax 2 



4x4 a 31 1 

— 4x log 



2xlog 
- 2x 



1 - Ax 2 

Ax 6 
I -Ax 2 



1 - 2x 
l + 2x ~ 

' -3.x 2 



1 



— 6x 5 log 



1 + 2x 12 
1 _ 9r 1 r 2 

- fl r 

1 + 2x 80 4 



+ 



(2.27) 



and 



h(x) = AI A (x) + BI B {x) + CI c (x) + 



(2.28) 



Here 



lA(x) = M r dy E^n r 

J~a It 2 JO 



dz X 



1 + 



V 
C(x,y) 



z(z + y) -C(x,y) 
E(x,y) r-y 



x x + y~\ 2x 2z — 1 
2n t + y - + 

a — z o + z ) y 1 + y 



, (2.29) 



lB(x) . £ 



dz 



C(x,y) 



and 





r * 




n < 













1 + z(z + y) -C(x,y) 
2 \ 2y + z + a 



x 



-i(2x + y)( a ; 2 + (a; + i/) 2 ) +j^(y + 2z) ^y + z + ^j + ( z ~ \) ' 



(2.30) 



y 2 jo 



dz 



1 + 



C(x,y) 



n 



y z — a 

2y + z + a 



z(z + y) - C(x,y) 

2 / 1 x 4\ 



X 



- — - \x A + x 2 (z - + (z - - 



lY 



■ + b \\ y + z + l) + { y + z + 2 l ) ( x + y) 2 + ( x + y) 4 )j 

--{2x + y) (x A + x 2 {x + yf + (x + yf) 

■jl-to + 2.) ((„ + * + i) 4 + (» + * + (* - i) 2 + (« - i)' 



(2.31) 



where the functions and C are defined as 



E(x,y) 



b{a + y) 



and 



C(x,y) 



M 2 6(a + y) - (m 2 + a)(l + y) ' 



(m 2 + a) (a + 
M 2 6(a + y) - (m 2 + a)(l + y) ' 



(2.32) 



(2.33) 



We may fit the integrands of I A {x), Ib{ x )i an d Ic{ x ) empirically after performing the z- 
integral: 

Integrand of Ia,b,c( x ) = fiA,B,c(v> r ) + f ) + fsA,B,c(y,r) (x + , (2.34) 



We return to Iq(x) and h(x) given by Eqs. (|2.14|) and ( |2.15|) , respectively, that may be written 
in the form where r = M 2 /(m 2 + a) and details of the functions fx and f% are given in the 
Appendix. This procedure produces equations that may be solved sequentially at each order 

— — — 2 

of x for the unknown parameters A, B, C, ■ ■ -, r, and g = m 2 +a , 



0(x°) 



r - 4 



3tt 



A A B C 

-4 + A-\ 1 

12 80 



+ 



9 2 



AF A0 + BF B0 + CF ( 



0{x 2 



Qtt 2 

(r - A) A - 4r 



CO 



(2.35) 



'3tt 



-16A + 



W C 



+ 



9 2 



3 5 

A(F A2 - AF A0 ) + B(F B2 - AF m ) + C(F C2 - AF C0 ) 



0(x 4 



6tt 2 

(r - 4)5 - 4r^ 



(2.36) 



3tt 



64A 



325 + 4C] 



+ 



67T 2 



A(F A4 - AF A2 ) + S(F B4 - AF B2 ) + C7(F C4 - 4F C2 ) 



(2.37) 



0(o; 6 ) 



(r - 4)C - 4rfi = - 



3tt 



256A 1285 - 

1 48C 

15 3 



+ 



6n 2 



(2.38) 



A(F M - 4F A4 ) + £(F B6 - 4F m ) + C(F C6 - AF C4 ) 

To solve to order x 2 , we neglect the terms in the expansion of the wave function to higher 
order that involve B, C, etc. This yields two equations from ( |2.35| ) and ( |2.36| ) in terms of A, 
r, and g. The solution to order x A involves neglecting terms C, etc., in the wavefunction and 
using ( |2.35D -( [2~3"TD , and so on. 

We can solve Eqs. (p . 35|) - ( p . 3 8|) numerically for r in terms of g. We show the result in 
Fig.l. For small values of m/g, we shall show that in the present model M/g approximately 
satisfies Eq. ( |2.25|) . 



III. VALIDITY OF THE APPROXIMATION. 



In this section we consider the contributions of Ii(x) and l2{x) or alternatively Ia{%) 
O, I B (x) (lOoP , and J c (x) ( ggg ) to the ground state mass M. 

At first, we neglect I\{x) and ^O 3 ?) in Eq. ( |2.13| ). We can solve the above equations ( |2.35| )- 
( |2.88| ) numerically for r in terms of g when we neglect these contributions of the integrals 
Ia{x), Ib{x), and Ic(x), ■ ■ Multiplying Eq. ( |2.13| ) by 1 — Ax 2 and setting D = . . . = for 
simplicity, we have 



O(x ) 
0(x 2 ) 
0(x 4 ) 



r - 4 



3tt 



-A + A-\ 1 

12 80 



Aq 

(r - 4)A - 4r = - 

37T 



— SB C 
-16A + — + - 
3 5 



(r - A)B - ArA 



3tt 
4q 

(r - 4)C - 4r£ = 

3tt 



6AA 



32B + AC 



256A 1285 
15 + 3 



48C 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



Then fl3.1|) to (|3.2|) with C — 0, after eliminating r and g, yield 



12A 3 + 144A 2 + 256A - 44AS + A 2 B - 33QB - 5 2 = 0. 



(3.5) 



Solving the above equation for A, we obtain three roots Ai(B), 

A2(B), and A^{B). The solutions where A and I? are real is shown in Fig. 2. The 
requirement that M is real rules out the solution A 3 in Fig. 2. 

We can express M 2 and m 2 in terms of A, B with the aid of ( |3.1|) , ( |3.2|) and (|3.3|): 



M 2 
m 2 



A 25 5 \ / / 1 A 3B B 
~A~3A + A8)/ \ 2 + 16~ 1QA + 192 



(3.6) 



We may also repeat this calculation, by including fl3T4|) , to find C. 

Our range of m/g is limited to — 2/^/3n. We have included these results in Fig. 1. For 
small values of m/g, M/g approximately satisfies (|2.25|) (dashed, lowest line). 
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We may repeat the procedures given above, but including the integrals Ia(%), Ib(%), and 
Ic(x). The changes from the solution of Sec. 3 and those calculated immediately above are 
under 10% for all the cases. Overall, Fig. 1 shows that M/g vs. m/g differs at most by less 
than 20% from the values obtained from the simplest model. 

IV. COMPARISON TO OTHER PARAMETRIZATIONS. 

We compare the relation M/m when different wavefunctions are used for bo(x). We begin 
by writing 't Hooft's equation for a two particle state || in the present notation, 

ttM 2 , f 1 1 \ (urn 2 \ r\ b (y) 



Eq. ( |4.1| ) corresponds to our Eq. ( j2.13j ). Both sides of Eq. are integrated with respect to 



x, and we obtain 



M 2 I * b (x)dx = 4m 2 r dx , h °^ X \ x . (4.2) 

i (1 — Ax 2 ) 



i 

2 "2 



When we substitute 't Hooft's wavefunction sin(n7r(x + 1/2)) = cos (nirx + — 2 1 - >7r ) in (|2.28| ) 



we find by inspection that only the odd numbers of n contribute to Eq. (|4.2j ) and obtain 

2M 2 /nn = 2m 2 Si(nvr), n = 1, 3, 5, • • • , (4.3) 
where the sine integral Si(z) is 

Si(z)= ( Z dt—. (4.4) 
Jo t 

For n — 1 

M 2 = 7rSi(7r)m 2 = 5.82 m 2 , or M/g = 2Alm/g. (4.5) 
which is in good agreement with ( |2.25| ). 
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If a wavefunctions of the form |7j b (x) = (1 — Ax 2 ) 13 is used, where (3 > and for the range 
of small P, we obtain bo(x) ~ 1 — 4j3x 2 . Upon comparison with ( |2.20| ) we get A = —4/3. One 
has for such wavefunctions (1 — Ax 2 ) 13 the boundary condition [?J 



— t r - 1 + n3cot(7ip) = 0. 

In the limit N c — > oo, we get tan7r/3 = 7T/3, which holds when (3 is small. 
For small (3 and N c = 3, Eq. ( |4.6|) is given by 

37rm 2 /V - tt 2 /3 2 (1 + 37rm 2 /8^ 2 )/3 = 0, 



(4.6) 



(4.7) 



so that 



A = -4(3 = -6m(l + 37rm 2 /8g 2 )- 1/2 /y/irg. 



We substitute Eq. flE^) in Eq. flggg ) and obtain 



M 2 



(1 + Svrm 2 ^ 2 ) 1 / 2 - m/2^g 
[1 + 3Trm 2 /8g 2 y/ 2 _ Sm/Ay/ng 



with the aid of < m/g < 2/y / 37r. We calculate Eq. (|4.9| ) and find 



m M m 
2.45— < — < 2.56—. 
9 9 9 



(4.8) 



(4.9) 



(4.10) 



The average slope of the numerical result by Hornbostel et al. ||, after a correction, is 



M/g = 3.1 m/g. 



(4.11) 



The SU(2) lattice gauge result of Hamer|| is 



M/g = 3.2m/ g. 



(4.12) 



A comparison of Eqs. (|2.25|) , ( pLop , (|4.10|) , Q4.11D and (|4.12|) , indicates that the ratios differ by 
about 30%. 
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Fig. 3 shows the comparison of the results of these models. It is apparent that all the 
models give values that are in general agreement. Thus, our simplest model b (x) = 1 + Ax 2 
is enough to obtain the approximate lowest ground state mass. 

V. SUMMARY AND CONCLUSIONS 

We have studied the ratio of the ground state mass M to the quark mass m. The 
wavefunction b Q (x) of ( |2.13 ) is expanded in a power series of even powers of x. We then 



matched equal powers of x on both side of the equation (|2.13|) . We have examined how each 



term of ( |2.19| ) affects the final result. When bo(x) is set equal to 1, we find to order x 2 that 



M = m = 0. 

Next, when bo(x) = 1+Ax 2 , the x° and x 2 equations yield the result Q2.25| ), M/g = \/6m/g. 



The substitution of b (x) = smnii(x + 1/2) in 't Hooft's equation gives for the ground state 
(n=l), M/g = [7rSi(7r)] 1 / 2 m/(7 = 2.41 m/g, where Si(x) is the sine integral. 

The simplest model of the wavefunction seems adequate to describe the features of the two 
dimensional light front description of mesons. The additional terms of the wavefunction bo{x) 
provide small corrections. 
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APPENDIX. DESCRIPTION OF THE FIT FOR I A (x), I B {x) AND I c {x) 
In general, our functions fiA,B,c,- an d fzA,B,c, - are parametrized in the form 

^ r) = ? Wr) S^H=' « Al > 

with different values of the parameters as shown in Table 1. The r-dependence of bij t k(r) is 
quadratic. 



Table 1: Values of parameters 





n 


7 


5 


V 


q 


flA 


3 





1 


1 





flB 


8 


1 











fic 


10 


1 








5 


Ha 


8 


1 


1 


2 





Hb 


3 





1 


1 





he 


6 


1 








5 
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Figure Captions 

Fig. 1: The solutions to the integral equations for M/g vs. m/g in the present paper with 
the wavefunction expanded to order x 2 , x 4 , and x 6 . The dotted line indicates the relation 
M/g = y/Em/g. 

Fig. 2: To order x 4 , the wavefunction is bo(x) = 1 + Ax 2 + Bx 4 . When r and g are 
eliminated from the matches, the parameters A and B are related. The real solutions are 
constrained as shown. 

Fig. 3: Results of the present work are compared to others calculations. Filled dots 
indicate the results of the present work, circles show the Hornbostel et al.'s results ||, open 
squares present the results of 't Hooft's large- N expansion 0, and diamonds indicate Sugihara 
et aPs results @. The masses M and m are given in units of J g 2 N c /2n. 
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<>(♦): 1 +Ax 2 +Bx 4 +Cx 6 with(out) l A;B; c" 
□ (■): 1+Ax 2 +Bx 4 with(out) l AB 



o (•): 1+Ax 2 with(out) l A 



1 - 



3 



0-o 







0.2 0.4 

m/g 

Fig. 1 



0.6 



A 
Fig. 2 



